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We report on a recent calculation of the electroweak ^(a 3 a s ) corrections to 3-jet production 
and related event-shape observables at e + e~ colliders. The calculation properly accounts for the 
experimental photon isolation criteria and for the corrections to the total hadronic cross section. 
Corrections to the normalised event-shape distributions, which are exemplarily discussed here for 
the thrust distribution at LEP and linear-collider energies, turn out to be at the few-per-cent level 
and show remnants of the radiative return to the Z pole even after inclusion of appropriate cuts. 
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1. Introduction 

Precision QCD studies at e + e~ colliders rely on the measurement of the 3-jet production cross 
section and related event-shape observables. The deviation from simple 2-jet configurations is 
proportional to the strong coupling constant a s , so that by comparing the measured 3-jet rate and 
event shapes (see, e.g., Ref. [[I]]) with the theoretical predictions, one can determine a s . Including 
electroweak coupling factors, the leading-order (LO) contribution to this process is of order a 2 a s . 

Owing to recent calculational progress, the QCD predictions for event shapes [||, ||] and 3- 
jet production [Q, |5|] are now accurate to next-to-next-to-leading order (NNLO, a 2 a s 3 ) in QCD 
perturbation theory. Inclusion of these corrections results in an estimated residual uncertainty of 
the QCD prediction from missing higher orders at the level of well below 5% for the event-shape 
distributions, and around 1% for the 3-jet cross section. Using these results (combined |^] with the 
previously available resummed expressions), new determinations of a s from event-shape and jet- 
production data were performed, resulting in a considerable improvement of the theory uncertainty 
to 3% from event shapes [Q] and below 2% from jet rates A further improvement can be 
anticipated for the event shapes from the resummation of subleading logarithmic corrections 

At this level of theoretical precision, higher-order electroweak effects could be of comparable 
magnitude. Until recently, only partial calculations of electroweak corrections to 3-jet production 
and event shapes were available [[it]], which can not be compared with experimental data directly. 
In Ref. [ |TT| ] we have presented the first calculation of the next-to-leading order (NLO) electroweak 
(a 3 a s ) corrections to 3-jet observables in e + e~ collisions including the quark-antiquark-photon 
{qqy) final states. Note that the QCD corrections to these final states are of the same perturbative 
order as the genuine electroweak corrections to quark-antiquark-gluon (qqg) final states. Since 
photons produced in association with hadrons can never be fully isolated, both types of corrections 
have to be taken into account. In this short article we supplement the results of Ref. [|TI]], where 
the differential thrust distribution and the 3-jet rate are discussed at the Z-boson resonance, by 
showing results on the thrust distribution for higher LEP and linear-collider energies. An exten- 
sive discussion of event-shape observables and more details of our calculation will be published 
elsewhere. 



2. Corrections to jet observables 

Event-shape measurements at LEP usually rely on a standard set of six variables y, defined 



for example in Ref. [12]: thrust T, C-parameter, heavy jet mass p, wide and total jet broadenings 
Bw and Bj, and 2-to-3-jet transition parameter in the Durham algorithm F3. The experimentally 
measured event-shape distribution (do/dy)/(Jh a d is normalised to the total hadronic cross section 
Chad- I n the perturbative expansion, it turns out to be most appropriate to consider the expansion of 
this ratio, which reads to NNLO in QCD and NLO in the electroweak theory 

1 do _ / a s \dA / a s \ 2 dB / a s \ 3 dC /a \ d<5 7 / a s \ / a \ ddnw 



Chad dy \2n) dy + \2n) dy + \2n) dy + \2n) dy + V27T/ \2n) dy ' (2-1) 

where the fact is used that the perturbative expansion of Gh a d starts at order a 2 . The calculation 
of the QCD coefficients A, B, and C is described in Refs. pL H]. The LO purely electromagnetic 
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contribution Sy arises from tree-level qqy final states without a gluon. The NLO electroweak coef- 
ficient 5ew receives contributions from the G(OL) correction to the hadronic cross section and from 
the genuine &{a i a s ) contribution to the event-shape distribution, 



Chad = CO 

such that 



1 da /a s \dA /a s \/a\ dd^ 



Go dy \2tiJ dy \2nJ \2nJ dy 



+ E £ ?> (2-2) 



d5ew dd^ ^ dA 



7^-^,17- (2-3) 
dy dy dy 

yields the full NLO electroweak correction. Both terms are to be evaluated with the same event- 
selection cuts. As shown in the following, many of the numerically dominant contributions, espe- 
cially from initial-state radiation, cancel in this difference. 

In the experimental measurement of 3-jet observables at e + e~ centre-of-mass (CM) energy 
y/s, several cuts are applied to reduce the contributions from photonic radiation. In our calculation, 



we apply the criteria used in the ALEPH analysis [12]. Very similar criteria were also applied by 
the other LEP experiments. Particles (including b quarks) contribute to the final state only if they 
are within the detector acceptance, defined by the production angle relative to the beam direction, 
|cos 1 < 0.965. Events are accepted if the reconstructed invariant mass squared s' of the final-state 
particles is larger than s cut = 0.8 Is. To reduce the contribution from hard photon radiation, the 
final-state particles are clustered into jets using the Durham algorithm with resolution parameter 
Jcut = 0.002. If one of the resulting jets contains a photon carrying a fraction z.y > Zy cat = 0-9 of 
the jet energy, it is considered to be an isolated photon, and the event is discarded. The event-shape 
variables are then computed in the CM frame of the final-state momenta, which can be boosted 
relative to the e + e~ CM frame, if particles are outside the detector acceptance. 

In the computation of the G[a) corrections to the total hadronic cross section, we include 
the virtual electroweak corrections to qq final states, and the real radiation corrections from qqy 
final states, provided the above event-selection criteria are fulfilled. The corrections to the event- 
shape distributions receive contributions from the virtual electroweak corrections to the qqg final 
state, the virtual QCD corrections to the qqy final state, and from the real radiation qqgy final 
state. To separate the divergent real radiation contributions, we used both the dipole subtraction 



method [|13J, |14J] and phase-space slicing [|15p, resulting in two independent implementations. Soft 
singularities are regularized dimensionally or with infinitesimal photon and gluon masses; they 
cancel in the sum of virtual and real corrections. Collinear singularities from initial-state radiation 
(ISR) are only partially cancelled. The left-over collinear ISR singularity is regularized by the 
electron mass and absorbed into the initial-state radiator function, which we consider either at 



fixed order, or in a leading-logarithmic (LL) resummation Q16[]. Owing to the specific nature of the 
event selection, collinear divergences from final-state radiation (FSR) are only partially cancelled. 
The left-over FSR singularity arises from the isolated photon definition, which vetoes on photon 
jets with z.y > Zy, cut- This singularity is absorbed into the photon fragmentation function, which 
we apply in the fixed-order approach of Ref. [Fl7]]. For the non-perturbative contribution to this 



function, we use the &{a) two-parameter fit of ALEPH [18]. The fragmentation contribution 



derived in Ref. [|17|] is based on phase-space slicing and dimensional regularization. We recomputed 



this contribution using subtraction and mass regularization [14]. 
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The one-loop diagrams are generated with FeynArts [19]. Using two independent inhouse 
Mathematica routines, one of which builds upon FormCalc p0[], each diagram is expressed 
in terms of standard matrix elements and coefficients of tensor integrals. The tensor coefficients 



are numerically reduced to standard scalar integrals using the methods of Ref. Q22|]. The scalar 
master integrals are evaluated using the methods and results of Ref. [p3Q, where UV divergences 
are regularized dimensionally. For IR divergences two alternative regularizations are employed, 
one that is fully based on dimensional regularization with massless light fermions, gluons, and 
photons, and another that is based on infinitesimal photon and gluon masses and small fermion 



masses. The loop integrals are translated from one scheme to the other as described in Ref. [24] 



The Z-boson resonance is described in the complex-mass scheme y25|], and its mass is fixed 
from the complex pole. The electromagnetic couplings appearing in LO are parametrized in the 
scheme, i.e., they are fixed via a = (Xg^ = V^G^M^ (l — M^/M|) /n. As the leading electromag- 
netic corrections are related to the emission of real photons, we fix the electromagnetic coupling 
appearing in the relative corrections by a = 05 (0), which is the appropriate choice for the leading 
photonic corrections. Accordingly the cross section for e + e~ — > qqg is proportional to cc G ^a & while 
the electroweak corrections to this process are proportional to a(0)a G 0$. The precise numerical 
input of our evaluations can be found in Ref. [0]. 

We performed two independent calculations of all ingredients resulting in two independent 



Fortran codes, one of them being an extension of Pole [21]. 



3. Results on the normalized thrust distribution 

In Fig. |l|, we display the differential thrust distribution normalised to the total hadronic cross 
section at different CM energies including NLO electroweak corrections and the relative corrections 
separately. The distributions are weighted by (l — T), evaluated at each bin centre. The corrections 
are obtained according to the expansion ( |2.1| ), retaining only terms up to LO in a s . The Born 



contribution is given by the A-term of (2.1), while the full G[a) corrections contain the tree- 
level qqy contribution 8 r and the NLO electroweak contribution 5ew (2-3). With "weak 0(a)" 
we denote the electroweak NLO corrections without purely photonic corrections, and "h.o. LL" 
indicates the inclusion of the higher-order ISR effects. 



As discussed in Ref. Q11Q for the Z pole in more detail, large ISR corrections cancel upon 
normalizing the event-shape distribution to the hadronic cross section, resulting in electroweak 
corrections of a few per cent. Moreover, effects from ISR resummation are largely reduced as well. 
Note that the photonic corrections develop a distinctive peak structure of up to 9% in size inside 
the thrust distribution, an effect that is a remnant of the radiative return to the Z pole, which is 
suppressed, but not fully excluded, by the event-selection cuts. The purely weak corrections are 
below 0.5 per mille at the Z pole, and only grow to the per-cent level for the linear-collider energy 
of yfi = 500GeV. 

Data on event-shape distributions and jet cross sections have been corrected for photonic ra- 
diation effects modelled by standard LL parton-shower Monte Carlo programs. They can thus not 
be compared directly with the NLO electroweak corrections computed here. Incorporation of these 
corrections requires a more profound reanalysis of LEP data, in order to quantify the impact of 
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Figure 1: Differential thrust distribution normalised to Chad at different CM energies y/s. 



the NLO electroweak corrections on precision QCD studies, such as the precise extraction of the 
strong coupling constant at NNLO in QCD [|7p. 
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